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Abstract. We investigate the Peierls transition in the one-dimensional Peierls- 
Hubbard model at half filling in the adiabatic approximation for the lattice. 
Depending on the value of the electron-lattice coupling constant g the equilibrium 
dimerization can be either enhanced or suppressed by the Hubbard interaction 
U. Applying second order perturbation theory we determine the critical value 
g c = 0.689348 below which the Hubbard interaction enhances the dimerization. 



PACS numbers: 71.10.Fd,71.30.+h 

1. Introduction 

The coupling of a one-dimensional metal to an elastic lattice results in an instability 
towards a lattice distortion known as Peierls transition pQ. This phenomenon can be 
observed in various quasi-one-dimcnsional materials, e.g. conjugated polymers |2] like 
polyacetylene, charge-transfer salts [3], or spin- Peierls compounds like CuGeC>3 . 

From a theoretical point of view, a first step towards a quantitative description of 
the Peierls transition has been made by Su, Schrieffer and Heeger (SSH) T . However, 
in their model neither the dynamics of the lattice nor the electron-electron interaction 
has been taken into account. Subsequently there has been made some effort to include 
phonons beyond the adiabatic approximation. While in the approach of SSH the 
Peierls transition occurs for any nonzero value of the electron- phonon (el-ph) coupling, 
a minimum strength of the el-ph coupling is required for the transition if quantum 
phonons are included 

In this paper we consider the limit of a static lattice and focus on the role of the 
electron-electron interaction. Correlation effects on the Peierls transition have been 
studied using a variety of methods including variational wave- functions [H||7], Hartree 
Fock plus perturbation theory jS] , quantum Monte Carlo [5j , numerical diagonalization 
of small systems [101 1111 1121 113) . bosonization ^] E3, an d incremental expansion 
[TBI IT7| . It is now generally accepted that for small values of the el-ph coupling 
constant g < g c the dimerization is enhanced when the Hubbard interaction U is 
switched on while for <7 > <7 C it is suppressed. Estimates for the critical el-ph coupling 
that have been obtained in the literature are g c « 0.75 in [Jj, g c w 0.76 in ^2] and 
g c w 0.69 in |16l I17j . For g not too close to g c the dimerization reaches a maximum 
at U of the order of the bandwidth and goes to zero for U — > oo. 

Apart from this qualitative agreement there is still a lack of rigorous results, in 
particular in the limit of weak el-ph coupling g <C g c which is difficult to address in 
the aforementioned approaches. In order to remove this uncertainties we study the 
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Peierls-Hubbard model in the weak coupling regime using second order perturbation 
theory with respect to the Hubbard interaction. This allows use to derive rigorous 
results in the limit C/ — > for all values of the el-ph coupling. 

The paper is organized as follows: In Section two we introduce the Peierls- 
Hubbard model and discuss the validity of the perturbative approach. The main 
results concerning the Peierls transition are presented in Section three. In Section 
four we discuss our results and compare them with different approaches. 



2. Model and perturbation theory 

We consider the one-dimensional Peierls-Hubbard model 

H = - t}^(l - a(u i+ i - tti))(cj a .Ci+i,<7 + h.c.) 

i,a 

+ U22,ni\mi + — ^{ui -u l+ if (1) 

i i 

where c\ a (ci a ) creates (annihilates) an electron with spin cr(=til) on a lattice of N 
sites, flic = c ia Ci a , and Uj is the deviation of the i'th atom from its equilibrium position 
in units of the lattice constant a. U and K are the on-site Hubbard interaction and 
the elastic constant of the lattice, respectively, and t is the hopping parameter. In Q 
the phonons are treated in the adiabatic limit which is formally obtained by letting 
the mass of the atoms go to infinity. Since the lattice distortion is assumed small 
compared to the lattice constant, the modification of the hopping amplitude can be 
treated in linear order via the parameter a. In the following we restrict ourselves 
to half filling where due to 2k p = ir/a the Peierls instability leads to an alternating 
lattice distortion 

Ui = (-l)*u (2) 

It is convenient to express the el-ph coupling and the dimerization amplitude by the 
dimensionless parameters g and 6, respectively, which are defined by 

9 = a ^/^ ( 3 ) 
6 = 2au (4) 
In terms of these parameters Hamiltonian reads 

H = -t]T(l + (-ir<S)( C U +1 , CT + h.c.) +UJ2^U + ( 5 ) 

We consider t, U and g as given model parameters while the dimerization 5 is free 
to adjust itself such that the total energy is minimal. In the following we measure 
all energies in units of the hopping amplitude t (i.e. we set t = 1). The ground 
state energy of the Hamiltonian JSJ can only be calculated analytically in the limiting 
cases (5 = and 6 = 1, respectively. S = corresponds to the ordinary Hubbard 
model which has been solved using Bethe ansatz ^S] and for 6 = 1 the system is 
completely dimerized. For other values of 5 one has to resort to approximations or 
to numerical methods. Since we are only interested in the weak coupling regime we 
use second order perturbation theory with respect to the Hubbard interaction. First 
of all we want to justify the assertion that perturbation theory can indeed be applied 
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for the model J5J|. In the case 6 = the exact ground state energy known from the 
Bethe ansatz solution j^] can be expressed as an integral over Bessel functions. As 
pointed out by Economou and Poulopoulos [15] this integral can be converted into 
an asymptotic series in U. Metzner and Vollhardt [201 showed that the exact second 
order term of this asymptotic expansion can be calculated using perturbation theory. 
Furthermore, they conjectured that this is true not only for the second order term but 
for all coefficients of the expansion. For 6 > the conditions for applying perturbation 
theory are more favorable than for 6 = since the unperturbed ground state energy Eq 
is separated from all excited states by an energy gap A = 48 which avoids divergent 
contributions from small denominators. We expect therefore that the ground state 
energy of the Peierls-Hubbard model can be expanded in a series with a finite radius 
of convergence, Uq(5) for all 5 > and that the coefficients of this expansion can be 
obtained from perturbation theory. Krivnov and Ovchinnikov |21| have pointed out 
that within the parquet approximation for a continuum version of the model JSJ the 
leading n'th order contribution to the ground state energy is oc 8 2 U n \u n+1 8. Similar 
conclusions were obtained by Horovitz and Solyom |22| using a scaling approach and 
by Kivelson et al |23| by comparison of exact results on the massive Thirring model 
with the perturbation theory for a model of spinlcss fermions. Fermions with spin were 
considered in |24| and again a perturbation expansion in powers of U In 6 was obtained. 
Thus the radius of convergence behaves as Uq(5) oc l/\ In <5| and even for exponentially 
small dimcrization 6 perturbation theory can be applied over some finite range of U. 
Writing e(6, U) = E(6, U)/N, the ground state energy per site of the electronic part 
of the expansion of e in powers of U takes the form 

e(5, U)-j = e (6) + e 2 (6)U 2 + e 4 (6)U 4 + ... (6) 

where we have used that due to particle-hole symmetry e(6, U)—U/4 is an even function 
of U. In particular, in the case of complete dimerization, 6 = 1, it is straightforward 
to calculate the exact ground state energy and expand it in a series which converges 
for U < 8: 

U I U 2 U 2 C/ 4 

e ( 1 '^-4-- 2 V 1+ 64-- 2 - 64 + 16384 + --- " (7) 

3. Equilibrium dimerization 

In order to determine the equilibrium dimerization we minimize the ground state 
energy of the Hamiltonian (J3J), including the lattice contribution, with respect to 6 
and obtain the condition 
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e' Q (6) + e' 2 (6)U 2 + e^t/ 4 + ... + — = (8) 

9 

Eq. |(HJ defines the function 6{U) which can as well be expanded in a series containing 
only even powers of U 

6{U) = 6 + 6 2 U 2 + 6 4 U 4 + ... (9) 
Reinserting Eq. into Eq. JSJ and comparing order by order in U yields 

5 = -g 2 e' (6 ) (10) 
4(M 



(11) 



d 4 = „_ 2 1 ^ ( i2 ) 



g- 2 + ^(6 ) 

±6 2 e'>>( 6 )+ 6 2 e>>(6 ) + e' 4 (6 ) 

g- 2 + 4(6 ) 
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Figure 1. Second order contribution to the ground state energy per site, 



6 = eo + e 2 U 2 + EiU 4 
minimum occurs at S c 



- as function of the dimerization parameter 8. The 
: 0.310523. 



The zeroth order dimerization has been calculated in 0j . In the limit g <C 1 it is given 

by 

5o =4exp(-^-l) (13) 

Now we turn to the second order contribution oc U 2 of the ground state energy. 
Considering the Hubbard interaction in Eq. (jSJ) as perturbation, the second order 
coefficient of the series ^ reads 

1 ^ \{n\D\0)\ 2 



E n 



En 



(14) 



where E n is the energy of the eigenstate \n) of the noninteracting Hamiltonian Hq, 
and D = ^\ Uifny. Representing |n) in the one-particle basis of Hq, Eq. (|14|) can be 
converted into 

, , .... , dk' 

e 2 (<5) 



1 

32 



dk 
2^ 



2tt 



2tt 



h _|_ e i{<Pk+q—<Ph+<Pi.' 



e am 



where 



£fc = 2\/ cos 2 k + S 2 sin 2 fc 



G = 



1 + e 



-2iA: 



+ (5(1 - e" 2ife ) 



-(15) 
(16) 

\1 + c- 2lk + 5(1 ~ e~ 2tk )\ ' ^ 
Eq. (|15J) is equivalent with Eq. (C5) in Ref. [B] in the limit of zero staggered 
magnetization. The integral (|15fl can be calculated numerically with very high 
precision. The result shown in Fig. ^ agrees with the analytical values £2(0) = 
-(7/16tt 3 )C(3) = -0.016961 and e 2 (l) = -1/64 = -0.015625. More importantly, 
€2(6) has a minimum at S c — 0.310523 which means that according to Eq. (|ll|l the 
sign of the coefficient 62 changes from plus to minus. According to Eq. (fltlfl the critical 
dimerization S c corresponds to a critical el-ph coupling constant 



9c = 



K(Sc)\ 



= 0.689348 



(18) 




Figure 2. Second order coefficient S2 of the equilibrium dimerization S(U) = 
So + S2U 2 + S4U 4 + ... as function of the el-ph coupling constant g. Note that 82 
changes sign at g c = 0.689348 which is related to <5 C = 0.310523 via Eq. 1181 . 



Therefore for g < g c the dimerization is enhanced by the Hubbard interaction while 
for g > g c it is suppressed. The second order coefficient 82 of Eq. © which is displayed 
in Fig. |21 has a maximum at g ~ 0.464. In the limit of small el-ph coupling, g <C 1, 
Eq. QII \ simplifies to 

82 = ~e' 2 (S ) (19) 

From a careful analysis of the numerical data for small values of 5 we conjecture that 

e 2 (S)-e 2 (0)^S 2 (ln^y (20) 

with c w 1.5, in agreement with the result of the parquet summation |21| . Inserting 
this expression into Eq. I|19() and using the asymptotic formula for So(g) given in Eq. 
l|13fl we obtain the result 

£ oc „"» (21) 

i.e. the relative weight of the second order term in the expansion of the dimerization 
6(U) = 5q + 82U 2 + S4U 4 + . . . increases strongly for g — * 0, although both So and 62 
are exponentially small in this limit. 

Unfortunately it becomes rather complicated to calculate the fourth order term 
£4 (S) using perturbation theory. We have therefore chosen a different approach based 
on the exact numerical diagonalization of small systems (up to L — 14 lattice sites) 
using the Lanczos algorithm which allows us to calculate the ground state energy with 
a relative precision of better than 10~ 14 . In order to avoid even-odd oscillations we 
have implemented antiperiodic boundary conditions for L — An and periodic boundary 
conditions for L = An + 2. This guarantees that the ground state is non-degenerate. 
Since we know the exact second order term e 2 {8, L) (for given system size L) it is 
possible to extract €4(6, L) with very high precision. However, one is left with the 
problem of extrapolating these values to L = 00. We have tried various extrapolation 
procedures. For 6 > 0.2 the ansatz 64(6, L) — £4(5,00) + aexp(— bL + c/L), which 
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Figure 3. Fourth order contribution of the ground state energy per site, 
e(L) = £o(L) + f-2(L)U 2 + £4 (L)U 4 + . . . , as function of the dimerization parameter 
S for L = 10, 12, 14 and oo. The dashed curve is meant as a guide to the eye, only. 
The solid line is the asymptotic result t4(<5) = 2 _14 (1 + 3(1 — <5)/2) valid in the 
limit S — > 1. 



is motivated by the existence of a gap in the excitation spectrum, yields very stable 
results for £4 (5, 00). In particular, the asymptotic behavior close to the fully dimerized 
limit, £4(15) = 2~ 14 (1 + 3(1 — S)/2 + ...) for S — > 1 is correctly recovered (see Fig. 
|2J|. However, for 5 < 0.2 we could not perform a reliable finite-size scaling of our 
data. The reason is that the correlation length, which is inversely proportional to the 
gap, is getting much larger than the size of the systems that can be diagonalized 
numerically. On the other hand, for 5 = where there is no gap, an unbiased 
polynomial fit yields £4(0, 00) = —1.547 x 10~ 4 which is quite close to the exact 
value Q2I e 4 (Q) = -(93/2048tt 5 )C(5) = -1.538 x 1(T 4 . In Fig. we plot e 4 (0,i) for 
L = 10,12,14 together with the extrapolated value £4(0,00) (whenever available) as 
a function of 5. Comparing with the corresponding results for £2 (see Fig. ^| we see 
that £4 is typically two orders of magnitude smaller than £2, which indicates that the 
weak coupling approach should be quite accurate in the region U ~ 1. 

For g < g c the dimerization S(U) reaches a maximum at U = f/ max which is given 

by 



Umax — A / 



2S 4 



(22) 



when we neglect the terms of order ~ U e and higher in Eq. . Expanding L^ ax in 
powers of g c — g we obtain for g close to g c 

U m U9)-<9c-g) 112 (23) 

with 

K = ; 5 uf {5c \ Wl u « 8 - 7 (24) 

<7<4(<y(i + 9c^c)) 

to be compared with the value n = 8.25 obtained in Ref. using the incremental 
expansion. Note that the error in the determination of k is only due to the limited 
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Figure 4. Dimerization S as function of U for g = 0.4,0.45,0.5,0.55,0.6 (from 
bottom to top). The full lines (circles) represent the second (fourth) order 
approximation of 5(u) = Sq + &2U 2 + 84U 4 + .... The values of U where <5 
becomes negative indicate the breakdown of the weak coupling approximation. 



precision of 64. Nevertheless, the fair agreement with the result of indicates that 
the finite-size scaling of our Lanczos data yields reasonable results for 64. 

In Fig.0]the equilibrium dimerization 8 of the Peierls-Hubbard model is displayed 
as function of U for several values of the el-ph coupling g. The circles represent the 
dimerization including the fourth order term in U while the solid curves correspond 
to the second order approximation. For small values of g the maximum of 8 is more 
pronounced and occurs around U ~ 2. This is somewhat smaller than the value 
U ~ 3 obtained using the incremental expansion 16 . Furthermore the maximum of 
8 is expected to shift to larger values of U with decreasing g but of course the precise 
position of the maximum depends on the terms of order oc U 6 and higher which have 
been omitted in our weak coupling approach. 

4. Discussion 

We have investigated the Peierls Hubbard model in the weak-coupling regime using 
second order perturbation theory and exact diagonalization for the fourth order 
contribution. The application of perturbation theory is justified by the existence of 
a gap in the excitation spectrum which removes the problem of small denominators. 
On the other hand, it is also due to this gap that finite-size scaling of the Lanczos 
data becomes problematic for small dimerization 8 < 0.2 where the correlation length 
exceeds the size of the systems that can be diagonalized numerically. The central 
result of this paper is the determination of the exact value of the el-ph coupling 
g c = 0.689348 which separates the region where the Hubbard interaction enhances 
the dimerization (<? < g c ) from the region where the dimerization is reduced (g > g c ). 
The corresponding dimerization parameter is found to be 8 C = 0.310523. There have 
been many attempts to determine these parameters using different methods. Using 
the Gutzwiller variational wave-function Baeriswyl and Maki [7] estimated g c « 0.75. 
Exact diagonalization of small systems has been used to obtain 8 C ss 0.4 in JT] and 
g c « 0.76 in [T2\, respectively. Finally, the most precise values up to now have been 
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derived using the incremental expansion technique, which is based on the numerical 
diagonalization of finite chains. In E| 9c = 0.69 is obtained, which is very close 
to the exact value. 

As mentioned before, all approaches that rely on finite-size extrapolation of small 
systems suffer from severe problems in the limit of small dimerization, 5 < 0.1, which 
corresponds to values of the el-ph coupling g < 0.5. In contrast, the perturbation 
theory approach can be applied for all values of g. In particular, it is possible to 
determine the asymptotic behavior of €2(5) in the limit <5 — > including logarithmic 
corrections. From there it can be inferred that the relative importance of the second 
order term (compared to the leading term) in the expansion of the dimerization 
S(U) = So + S2U 2 + 54U 4 + . . . goes as ~ g~ 6 for g — > 0, although both So and S2 are 
exponentially small in this limit. On the other hand, considering absolute values, the 
correlation enhancement of the dimerization is less dramatic and essentially limited to 
the region 0.3 < g < g c with the maximum occurring at g = 0.464. The parameters 
most frequently used for polyacetylene lie in the range between g — 0.39 [7] and 
g = 0.57 4 which emphasizes the importance of correlation effects in real Peierls 
systems. 
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